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Heat  intercept  insulation  system 
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Heat  intercept  insulation  system 


min  power(n,I,x,T) 
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Heat  intercept  insulation  system 


min  power(n,  I ,  x,T) 

subject  to  n  G  { 1,2,...,  nmax},  /  G  Xn+1 

<Ti<  Ti+ 1,  i  =  1,  2, . . . ,  n 


n+1 

Xi  =  L,  >  0,  i  =  1,  2, . . . ,  n  +  1 

2  =  1 
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The  general  MVP  problem 


min 

xGX 


subject  to  C(x)  <  0, 
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The  general  MVP  problem 


min 

xGX 


subject  to  C(x)  <  0, 


where 

x  =  (xc,  xd)  e  3£nC  x 

=  Xc  x  Xd,  where 

c(xd )  =  {xc  G  :  i{xd)  <  A(xd)xc  <  u(xd)} 
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The  general  MVP  problem 


min 

xGX 


subject  to  C(x)  <  0, 


where 

x  =  (xc,  xd)  e  3£nC  x 
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Choices  other  than  pattern  search 


Some  methods  that  come  to  mind  are: 

MINLP  methods:  cannot  handle  categorical 
variables 
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Choices  other  than  pattern  search 


Some  methods  that  come  to  mind  are: 

MINLP  methods:  cannot  handle  categorical 
variables 

Search  heuristics:  huge  numbers  of  evaluations  and 
very  limited  convergence  theory 

■  Simulated  annealing 

■  Tabu  search 

■  Evolutionary  algorithms 

Other  methods:  SQP/direct  search  with  1 
categorical  variable 
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Generalized  pattern  searches 

INITIALIZATION  of  directions  and  step  size 
For  k  =  1,2,... 

SEARCH  a  finite  set  of  mesh  points 
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Generalized  pattern  searches 


INITIALIZATION  of  directions  and  step  size 
For  k  =  1,2,... 

SEARCH  a  finite  set  of  mesh  points 


End 


POLL  neighboring  mesh  points 


Pattern  Search  for  Mixed  Variable  Optimization  Problems  -  p.6/35 


Generalized  pattern  searches 


INITIALIZATION  of  directions  and  step  size 
For  k  =  1,2,... 

SEARCH  a  finite  set  of  mesh  points 

POLL  neighboring  mesh  pointsU . 

UPDATE  parameters: 


End 
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Generalized  pattern  searches 


INITIALIZATION  of  directions  and  step  size 
For  k  =  1,2,... 

SEARCH  a  finite  set  of  mesh  points 

POLL  neighboring  mesh  points 

UPDATE  parameters:  "*( 

■  Success:  Accept  new  iterate 

End 
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Generalized  pattern  searches 

INITIALIZATION  of  directions  and  step  size 


For  k  —  1,2,... 

SEARCH  a  finite  set  of  mesh  points 
POLL  neighboring  mesh  points 
UPDATE  parameters: 


Success:  Accept  new  iterate 
Lailure:  Refine  mesh 


End 
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Details  of  kth  POLL  step 

Mesh:  M \  =  {pk  +  A kDz  :  z  G  Z^}, 
Poll  scti  Pfc  =  \jPk  H-  ^kd  •  d  G  Dfc  C 
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Details  of  kth  POLL  step 

Mesh:  M \  =  {pk  +  A^Dz  :  z  G  Z^}, 
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where 

Pk  is  the  current  poll  center 
Ak  >  0  is  the  mesh  size  parameter 
-Dyt,  D  are  positive  spanning  sets 
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Details  of  kth  POLL  step 

Mesh:  M \  =  {p/,:  +  A^Dz  :  z  G  -2+^}, 

Poll  scti  Pfc  =  \jPk  H-  ^kd  •  d  £  Dfc  Cl  £}. 

where 

■  pk  is  the  current  poll  center 
Ak  >  0  is  the  mesh  size  parameter 
-Dyt,  D  are  positive  spanning  sets 

Examples:  D  =  [/,  — /]  D  =  [/,  — e 

One  of  these  directions 
should  be  a  descent  direction. 

Derivative  information  can 
reduce  poll  set  to  a  singleton 
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Details  of  kth  POLL  step 


Mesh:  M \  =  {pk  +  A^ Dz  :  z  E  }, 

Poll  scti  Pfo  =  \jPk  H-  ^kd  •  d  £  Dj$  C 

where 

Pk  is  the  current  poll  center 
Afc  >  0  is  the  mesh  size  parameter 
Dk ,  D  are  positive  spanning  sets 

Examples:  D  =  [/,  —  /]  D  =  [I,  — e 

One  of  these  directions 
should  be  a  descent  direction. 

Derivative  information  can 
reduce  poll  set  to  a  singleton 
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Definition  of  local  optimality 

x  =  (xc,  xd )  G  X  is  a  /oca/  minimizer  of  f  with  respect  to 
neighbors  N(x)  C  X  if  Be  >  0  such  that  f(x)  <  f(v ) 

V»£Xn  (J  (B(y\e)xyd). 

y£j\f(x) 
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Definition  of  local  optimality 


x  =  (xc,  xd )  G  X  is  a  /oca/  minimizer  of  f  with  respect  to 
neighbors  J\f{x)  C  X  if  Be  >  0  such  that  f(x)  <  f(v ) 

v«  e  xn  U  (B(f,e)xyd). 

y&M(x) 

Xc  x 


Xc  x 
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Heatshield  discrete  neighbors 

Replace  any  single  insulator  with  a  different  type 
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Heatshield  discrete  neighbors 


Replace  any  single  insulator  with  a  different  type 

Remove  any  intercept  with  its  left  insulator,  and 
increase  the  thickness  of  its  right  insulator  to  absorb 
the  deficit 

Add  an  intercept  at  any  position: 

■  The  existing  insulator  is  divided  (rounded  to  the 
mesh) 

■  The  cooling  temperature  is  set  to  the  average  of  the 
two  intercepts  adjacent  to  it  (rounded  to  the  mesh) 
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?r  A, 


Construction  of  the  poll  set 
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?r  ft- 


Construction  of  the  poll  set 


Pk  =  {u,v,  w} 
Af{xk)  =  {xk,y^,y^} 


j/l  G  A f(xk)  satisfies 

f{xk)  <  f(Vi)  <  f(xk)+€ 
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?r  ft- 


Construction  of  the  poll  set 


Pk  =  {u,v,  w} 

Af{xk)  =  {xk,y^,y^} 

Xk  =  {y\}  U  {a,  b,  c} 

y\  E  M{xk)  satisfies 

f(xk)  <  f{Vi)  <  f{xk)+€ 

Poll  Set:  Pk  U  A f(xk)  U  Xk 
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Filter  GPS  for  nonlinear  constraints 


h(x)  =  ||C'(x)+||2 
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Filter  GPS  for  nonlinear  constraints 


h(x)  =  ||CY(:c)+||2 


Poll  center  is  either  best 
feasible  point  or  least 
infeasible  point. 


Pattern  Search  for  Mixed  Variable  Optimization  Problems  -  p.  13/35 


Filter  GPS  for  nonlinear  constraints 


h{x)  =  ||C'(x)+||2 


Poll  center  is  either  best 
feasible  point  or  least 
infeasible  point. 

For  each  trial  point  x, 
h(x)  and  f(x)  are  plot¬ 
ted  on  the  bi-loss  map. 
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Filter  GPS  for  nonlinear  constraints 


h{x)  =  ||C'(x)+||2 


Poll  center  is  either  best 
feasible  point  or  least 
infeasible  point. 

For  each  trial  point  x, 
h(x)  and  f(x)  are  plot¬ 
ted  on  the  bi-loss  map. 

If  x  is  unfiltered,  it  is 
added  to  the  filter;  oth¬ 
erwise,  the  mesh  is  re¬ 
fined. 
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Construction  of  the  poll  set 


Pk  =  {u,v,  w} 

N{Pk)  =  {Pk ,  Vi ,  2/2 } 
=  {y{}  U  {a,  b,  c} 


Di  G  N{pk)  satisfies 
the  extended  poll  criteria 


Poll  Set:  Pk  U  Af(pk)  U  Xk 
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Local  filter  for  extended  polling 


Main  Filter 
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Local  filter  for  extended  polling 


Main  Filter  Local  Filter 
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Filter  GPS  algorithm  for  MVP 

INITIALIZATION:  Set  Ao,  £  >  0,  and  populate  filter 
For  k  =  1,2,...,  do 

Update  poll  center  pk  and  extended  poll  triggers  ^ 
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Filter  GPS  algorithm  for  MVP 
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For  k  =  1,2,...,  do 

Update  poll  center  pk  and  extended  poll  triggers  ^ 
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While  trial  points  are  filtered  do: 
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Filter  GPS  algorithm  for  MVP 

INITIALIZATION:  Set  Ao,  £  >  0,  and  populate  filter 
For  k  =  1,2,...,  do 

Update  poll  center  pk  and  extended  poll  triggers  ^ 

Compute  incumbent  values  f['\  //,  //( 

While  trial  points  are  filtered  do: 

SEARCH:  Any  finite  strategy  on  the  mesh 
POLL:  Evaluate  points  in  Pk  U  M{pk) 

■  EXTENDED  POLL:  Evaluate  points  in  Xk(£fk,  $) 
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Filter  GPS  algorithm  for  MVP 
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Update  poll  center  pk  and  extended  poll  triggers  ^ 
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While  trial  points  are  filtered  do: 
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■  EXTENDED  POLL:  Evaluate  points  in  Xk(£fk,  $) 
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Filter  GPS  algorithm  for  MVP 


INITIALIZATION:  Set  A0,  £  >  0,  and  populate  filter 
For  k  =  1,2,...,  do 

Update  poll  center  pk  and  extended  poll  triggers 

Compute  incumbent  values  fk\  h\. 

While  trial  points  are  filtered  do: 

SEARCH:  Any  finite  strategy  on  the  mesh 

■  POLL:  Evaluate  points  in  Pk  U  M{pk) 

■  EXTENDED  POLL:  Evaluate  points  in  Xk(£fk,  $) 

■  Update: 

If  (found),  set  Afc+1  >  A k  and  update  filter 
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Filter  GPS  algorithm  for  MVP 


INITIALIZATION:  Set  A0,  £  >  0,  and  populate  filter 
For  k  =  1,2,...,  do 

Update  poll  center  pk  and  extended  poll  triggers 

Compute  incumbent  values  fk\  h\. 

While  trial  points  are  filtered  do: 

SEARCH:  Any  finite  strategy  on  the  mesh 

■  POLL:  Evaluate  points  in  Pk  U  M{pk) 

■  EXTENDED  POLL:  Evaluate  points  in  Xk(£fk,  $) 


Update: 


If  (found),  set  Afe+1 
If  (not  found),  set  A 


>  Ak  and  update  filter 


k+ 1 


<  A 


k 
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Convergence  theory  assumptions 

All  iterates  lie  in  a  compact  set 
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Convergence  theory  assumptions 
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Convergence  theory  assumptions 


All  iterates  lie  in  a  compact  set 

The  linear  constraint  matrix  A  is  rational 

The  mesh  directions  conform  to  the  geometry  of  Xc 

All  discrete  neighbors  lie  on  the  current  mesh 

The  set- valued  neighborhood  function  J\f  :  X  — »  2X 
satisfies  a  notion  of  continuity. 
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Limit  points  of  the  algorithm 


3  subsequence  K  such  that  lim  =  0,  with  limit  points: 

keK 

1-  P  =  )im  pk,  where  pk  G 

2.  y  —  lim  ijk ,  where  Vk  e  Af(pk)  and  y  e  A f{p). 

keK 

3.  z  =  lim  Zyt,  where  Zfc  are  extended  poll  endpoints. 
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Limit  points  of  the  algorithm 


3  subsequence  K  such  that  lim  =  0,  with  limit  points: 

keK 

1-  P  =  )im  pk,  where  pk  G 

2.  y  —  lim  ijk ,  where  Vk  e  Af(pk)  and  y  e  A f{p) 

keK 

3.  z  =  lim  Zyt ,  where  Zfc  are  extended  poll  endpoints 


Filter  convergence  results 


Let  D(p)  be  the  set  of  polling  directions  used  i.o. 
h  continuous*  at  p  and  y  =>•  h(p)  <  Kv) 
f  continuous*  at  p  and  y  and  pk  =  pf  i.  o. 

f(p)  <  f(y ) 

h  Lipschitz*  nearp  =4>  h°(p;  (d,  0))  >  0  \/d  E  D(p) 

f  Lipschitz*  near  p  and  p/,  =  p^  i.  o. 

=>  /°(p;  (d,  0))  >  0  Vd  G  D(p) 

h  strictly  differentiable*  atp  and  =>  V//(p)  =  0 
Similar  results  hold  for  certain  5 
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Filter  convergence  results 


/  strictly  differentiable*  at  p  and  pk  =  pf.  i.  o. 

=4>  — Vc/(j3)  €E  C%.  Similar  results  hold  for  certain  z 


*  with  respect  to  the  continuous  variables 
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Filter  convergence  results 
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Heat  intercept  insulation  system 
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Heat  intercept  insulation  system 


min  power(n,I,x,T) 
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Heat  intercept  insulation  system 


min  power(n,  I ,  x,T) 

subject  to  n  G  { 1,2,...,  nmax},  /  G  Xn+1 

<Ti<  Ti+ 1,  i  =  1,  2, . . . ,  n 


n+1 

Xi  =  L,  >  0,  i  =  1,  2, . . . ,  n  +  1 

2  =  1 
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Previous  heat  shield  studies 


Hilal  &  Boom:  1-3  intercepts,  single  insulator  type, 
constant  cross-sectional  areas 
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Previous  heat  shield  studies 


Hilal  &  Boom:  1-3  intercepts,  single  insulator  type, 
constant  cross-sectional  areas 

Hilal  &  Eyssa:  1  -3  intercepts,  single  insulator  type, 
variable  cross-sectional  areas 

Kokkolaras,  Audet  &  Dennis:  Variable  number  of 
intercepts,  multiple  insulator  types,  constant 
cross-sectional  areas 

Current  work:  Variable  number  of  intercepts, 
multiple  insulator  types,  variable  cross-sectional  areas 
load-bearing  nonlinear  constraints 
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Nomenclature 


KT-h) 

Pi 

Ci 

Qi 

A 

e(T ;  I,) 

P{h) 

F 

^max 


=  Thermal  conductivity  function  for  insulator  i 
=  Power  applied  to  intercept  i 
=  Thermodynamic  cycle  coefficient  at  intercept  i 
=  Heat  flow  from  intercept  i  to  i  —  1 
=  Cross-sectional  area  of  insulator  i 
=  maximum  allowable  stress  function 
=  unit  thermal  expansion  function 
=  density  of  the  insulator  i  material 
=  load  (force)  to  be  placed  on  the  system 
=  maximum  allowable  mass  of  the  insulators 
=  maximum  allowable  %  thermal  contraction 
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Heat  shield  objective  and  nonlinear 
constraints 


n 


n 


Minimize  power:  Vj  Pi  —  YJ  Ci 


T 


H 


i—  1 


By  Fourier’s  law:  qi  = 


A 


2=1 

•T. 


-  1  )  (Qi 


x 


k(T ;  Ii)dT,  i  =  1 


i  JT 


i—  1 


9*-i) 

2? . . .  ?  Tt  H-  1 
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Heat  shield  objective  and  nonlinear 
constraints 


n 


n 


Minimize  power:  Vj  Pi  —  YJ  Ci 


T 


H 


i—  1  i=  1 

‘T, 


^7  -  1 )  (ft 


By  Fourier’s  law:  qt  = 


A 


a; 


fc(T;  Ii)dT ,  i  =  : 


i  —  1 


Stress: 


F 

A 


<  min{cr(T;Ii)  :  TV  1  <  T  <  V} 


-  ?*-i) 

L ,  2, . . . ,  77,  T- 1 

(binding) 
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Heat  shield  objective  and  nonlinear 
constraints 


n 


n 


Minimize  power:  Vj  Pi  —  YJ  Ci 


T 


H 


i= 1 


By  Fourier’s  law:  qt  = 


A 


i—  1 

'Ti 


X 


k(T ;  Ii)dT ,  i  =  : 


i  —  1 


Stress: 


F 

A 


<  min{(j(T;/j)  :  T;_i  <  T  <  Ti} 


n+1 

Mass:  '^^p(Ii)Aixi  <  m 

max 

i=l 


~  Qi- 1) 

l  ,  2, . . .  ?  a  i 

(binding) 
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Heat  shield  objective  and  nonlinear 
constraints 


n 


n 


Minimize  power:  Vj  Pi  —  YJ  Ci 


T 


H 


i—  1  i—  1 

'Ti 


=7  -  l )  (ft  -  ft-i) 


By  Fourier’s  law:  = 


A 


x 


k(T;Ii)dT ,  *  =  1,  2, . . . ,  n  +  1 


i  —  1 


Stress: 


F 

A 

n+1 


<  min{(j(T;  /j)  :  T;_i  <  T  <  T;}  (binding) 


Mass:  ^  z  ^  //7  rnax 


i—  1 


Contraction: 


tr  \  fc(r;  U)dT  ) 


<5 

-  100 
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Materials: 

Nylon 

Teflon 

304  Stainless  Steel 
1 020  Low  Carbon  Steel 


6063-T5  Aluminum 
Fiberglass  Epoxy  (normal  cloth) 

Fiberglass  Epoxy  (plane  cloth) 
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Materials: 

Nylon  6063-T5  Aluminum 

Teflon  Fiberglass  Epoxy  (normal  cloth) 

304  Stainless  Steel  Fiberglass  Epoxy  (plane  cloth) 

1 020  Low  Carbon  Steel 

Material  Data  from  Lookup  Tables  or  Graphs: 

Thermal  Conductivity  Unit  Thermal  Contraction 

Maximum  Allowable  Stress  (Tensile  Yield  Strength) 
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Heat  shield  implementation 

Materials: 

Nylon  6063-T5  Aluminum 

Teflon  Fiberglass  Epoxy  (normal  cloth) 

304  Stainless  Steel  Fiberglass  Epoxy  (plane  cloth) 

1 020  Low  Carbon  Steel 

Material  Data  from  Lookup  Tables  or  Graphs: 

Thermal  Conductivity  Unit  Thermal  Contraction 

Maximum  Allowable  Stress  (Tensile  Yield  Strength) 

Interpolation/Integration:  Cubic  splines,  Simpson’s  rule 

Search/Poll:  No  Search,  Poll  around  pf 
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Source 

Insulators 

Hilal  &  Boom 

68.6 

Ep  Ep  Ep 

Kokkolaras  et  al. 

25.3 

NNNNNNNEEET 

Kokkolaras  rerun 

25.59 

NNNNNNTEETT 

Hilal  &  Eyssa 

53.2 

Ep  Ep  Ep 

with  stress  constraint 

24.55 

L  '  L  1  Lj  1  L  '  L  1  L  1  L  1  Lj  '  L  '  Lj  1  L  1 

with  all  constraints 

23.77 

L  1  L  1  L  1  L  1  7?  Lj  1  L  1  L  1  L  1  L  1  L  1 
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shield  computational  results 


Source 

Insulators 

Hilal  &  Boom 

68.6 

Ep  Ep  Ep 

Kokkolaras  et  al. 

25.3 

NNNNNNNEEET 

Kokkolaras  rerun 

25.59 

NNNNNNTEETT 

Hilal  &  Eyssa 

53.2 

Ep  Ep  Ep 

with  stress  constraint 

24.55 

L  '  L  1  TP  L  '  L  1  L 1  L '  ZT1  L  '  L 1  L  1 

with  all  constraints 

23.77 

L  1  L  1  L  1  L  1  Lj  '  Z71  77*  L 1  L  1  Z7'  L  1 

Parameters: 

T h  —  300  K,  T a  —  4.2  K,  L  —  100  cm,  nmax  —  10 
F  =  250  kN,  mmax  =  10  kg,  5  =  5% 
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shield  computational  results 


Source 

PL  TVE1 

A  L  cm  i 

Insulators 

Hilal  &  Boom 

68.6 

Ep  Ep  Ep 

Kokkolaras  et  al. 

25.3 

NNNNNNNEEET 

Kokkolaras  rerun 

25.59 

NNNNNNTEETT 

Hilal  &  Eyssa 

53.2 

Ep  Ep  Ep 

with  stress  constraint 

24.55 

L  1  L  1  Z7»  L  '  77'  L  1  L  '  L  '  L  '  Lj  1  L  1 

JLJLJLJLJLJLJLJLiLiL-Uj 

with  all  constraints 

23.77 

L  1  L  1  L  1  L  1  Lj  1  L  1  L  1  L  1  L  1  L  1  L  ’ 

JLJLJLJLJLJLJLJLiLiL-Uj 

Parameters: 

T h  —  300  K,  T a  —  4.2  K,  L  —  100  cm,  nmax  —  10 
F  =  250  kN,  mmax  =  10  kg,  5  =  5% 

Termination:  A*,  <  .15625 
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Power  Required  Power  ReqUired 

ro  ro  fn  ' 


Profile  of  heat  shield  run 
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shield  filter  pr< 


Filter  after  150  evaluations 


Zoom  of  filter  on  left 


0  0.05  0.1  0.15  0.2  0.25  0.3 


Zoom  of  filter  on  left 
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Conclusions 


We  have  found  an  effective  method  for  mixed  variable 
problems...  BUT... 
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